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Abstract. We provide algebraic criteria for the unitarity of linear quan- 
tum cellular automata, i.e. one dimensional quantum cellular automata. 
We derive these both by direct combinatorial arguments, and by adding 
constraints into the model which do not change the quantum cellular 
automata's computational power. The configurations we consider have 
finite but unbounded size. 

1 Motivat ions 

One could say that the central question in theoretical computer science is 'What 
are the resources necessary for comptitation, or Information processing?'. Ulti- 
mately this question is dictated by the physical laws which surround us. Quan- 
tum computation science has risen from this bàsic idea. It considers computers 
as physical, and hence possibly quantum systems. At the theoretical level it was 
demonstrated for instance that polynomial-time integer factorization is possible 
with such systems, as well as the search for an element in a unordered list of size 
n in time &{^/n). 

Cellular automata (CA) are arrays of celis, each of which may take one in a fi- 
nite number of possible state. These evolve in discrete time steps according to a 
global evolution A - which itself arises from the application of a local transition 
function S, synchronously and homogeneously across space. A popular example 
is Conway's 'Game of Life', a two dimensional CA which has been proven to be 
universal for computation. 

It is clear that CA are themselves physics-like models of computations, as they 
describe a world of small systems interacting locally, according to translation- 
invariant laws. Therefore it seems natural to study their quantum extensions. 
Moreover no classical control is required in such models, since computation arises 
as an emergent behaviour of the quantum celis' interaction. This is a key advan- 
tage to have as it reduces the need for environment interaction, and hence may 
reduce decoherence, the principal obstacle to realizing a quantum computer. For 
these two reasons Feynman in his seminal paper about quantum computa- 
tion, has argued that the study of quantum CA may prové an important path 
to a realistic physical implementation of quantum computers. Another series of 
legitimate aims is to endow quantum computation with spatio-temporal notions. 



or even to provide a bridge through which computcr scicncc notions, such as uni- 
versality, may contributc to modern theoretical physics. To put it difïerently such 
works are a contribution to the understanding of dynamics in discrete, quantum 
spacetime - but froni an idcalized, computer-science viewpoint. 
Similarly to thcir classical counterparts linear quantum cellular automata (LQCA) 
consist of a row of identical, finita dimensional, quantum systems. These evolve 
in discrete time steps according to a global evolution A - which itself arises 
from the application of a local transition function 5, homogeneously and syn- 
chronously across space. But in order to grant LQCA the status of physically 
acceptable model of computation, one must ensure that the global evolution A 
is physically acceptable in a quantum theoretical setting, i.e. one must ensure 
that A is unitary. Unfortunately this global property is rather non-trivially re- 
lated to the description of the local transition function 6 ~ witness of this the 
abundant literature on reversible cellular automata (RCA), tackling the classical 
counterpart of this issue. It is actually a very surprising fact that so much has 
been done to study RCA - when reversibility is not so much of a crucial feature 
to have in classical computation. A frequently encountered argument states that 
all consumption-less, zero-heat micro-mechanical device need be reversible. But 
tracing back the origins of this argument, we find quantum physical considera- 
tions once more |ü]. 

One way to approach this issue is to find a decision procedure, which given S telis 
whether A is unitary. This test should be performed efficiently, so as not to carry 
any of the complexity of the computation, and be applied to any candidate S as 
a mean to exclude the non-physical ones. Such a strong contribution was indeed 
achieved by Dürr et al. It does not put an ending point to the problem however, 
because the number of local transition functions which do indeed induce a uni- 
tary global evolution is likcly to be rather scarce, as was the case for RCA 
Moreover this rclatively complicated decision procedure takes an elegant detour 
via finite automata, but one which does not guide us to understand which (5's 
will eventually yield a physical A. 

Physicists are used to checking whether an evolution is physically acceptable, 
but they like to do so algebraically (e.g. WU = UW — I for unitarity, when U 
is a finite matrix) . Much in the same way computer scientists are used to check- 
ing whether a program is vàlid, but like these criteria to be syntactic. When 
this is not the case, we tend to consider that the definition chosen to formal- 
ize the model of computation is in fact too lose. Indeed once universality has 
been reached, adding more cxpressiveness does not mean adding more conipu- 
tational power, but only more ways of expressing the same computation. An 
undesirable excess arises when the syntax proposed by the definition allows the 
description of non-valid (unrcalistic, non-physical) programs, thereby requiring 
that the user performs non-trivial (non-syntactic) decision procedures to exclude 
those instances. This is the current state of affairs with LQCA. 
Therefore a different, complementary approach is to tighten the definition of 
linear quantum cellular automata, i.e. to seek for a more restrictive definition 
whose unitarity may be checked algebraically/syntactically, and yet capable of 



expressing the exact same set of global evolutions as our original definition. In 
this paper we provide algebraic characterizations of unitary linear quantum cel- 
lular automata. We derive these both by dircct combinatorial arguments, and 
by adding constraints into the model which do not change thc quantum cellular 
automata's computational power. 

The breakdown of the core of this paper will be given after LQCA are presented 
formally, in the foUowing section. Our main theorem is stated in the conclusion 
section, and discussed in comparison to some related approaches. 
Notations. Throughout the paper we will denote by TLs the hilbert space whose 
canonical orthonormal base vectors are identified with the elements of the count- 
able set S. E.g. 'H{aa,ab.ba,bb} ís the four dimensional space with canonical or- 
thonormal base {|aa), \ah), \ba), \hh)}. This means that any vector a\aa) + I3\ab) + 
7|òa) + 5\bh) with a, /3, 7, (5 € C belongs to 'H{aa,ab,ba,bb}- Such a vector must be 
thought of as a superposition of the words aa, ab, ba, bb. Moreover the symbol 
is to denote the nuU vector, not to be confused with O the matrix containing 
only ones. 

2 The model 

We start with the definition proposed by jJJ and ^ . This definition will evolve 
throughout the paper. 

Working definition 1 (LQCA) 

A linear quantum cellular automaton (LQCA) is a 4-tuple A = {S,q, N,S), 
where (with qS = {q\ U S): 

- S is a finite set of symbols (i.e. "the alphabet", giving the possible bàsic states 
each celi may take); 

- q is a symbol such that q ^ S (i.e. "the quiescent symbol", which may be 
thought as a special state for empty celis ); 

- N is a set of n successive signed integers (i.e. "the neighbourhood" , telling 
which celi is next to whom); 

- S : 'H(qs)" ~^ Ti-qE is a function from superpositions of n .symbols words to su- 
perpositions of one symbol words (i.e. "the local transition function" , describing 
the way a celi interacts with its neghbours). 

Moreover S must verify: 

- the quiescent stability condition: [(5|g")) = \q)\; 

- the no-nullity condition: Vw G (qZ")", \5\w) 7^ O] . 

By 'successive' we mean that the number follow each other in unit step, i.e. the 
neighbourhood is an interval. In the literature these are sometimes referred to as 
simple neighbourhoods, but it is trivial to simulate non-simple neighbourhoods 
automata with simple neighbourhoods automata. At this point we need not have 
a normalization condition such as Vw € (qi^)", [||(5|iü)|| = l] . Configurations 
hold the bàsic states of an entire row of celis. As we will now formalize ours are 
finite but unbounded. Note that fixed-sized periòdic configurations^ as well 



as infinite configurations^^I have also been studied, leading to very difFerent 
results and proof mcthods (see Section|ülfor a discussion). 

Definition 1 (finite configurations, interval domains) 

A (finite) configuration c of the quantum cellular automaton A = {S,q, N,S) is 

a function c : Z *■ qS, with i i — > c{i) = Ci, such that there exists a (possibly 

empty) interval I verifying ï G / Cj G qS and i ^ I ^ Ci = q. Moreover we 
denote by idom{c) the smallest such interval I , referred to as the interval do- 
main of c. Then the interval [k + Tahi(N), í + niax(iV)] is denoted extidom{c) and 
referred to as the extended interval domain of c. For the all-quiescent configura- 
tion c = . . . qqq . . . we have idom{c) = 0, and we let extidom{c) = also. Finally 
the set of all finite configurations is denoted Cf, whilst the set of configurations 
having interval domains comprised within an interval J is denoted Cj . 

Definition 2 (indexing conventions) 

Given a configuration c of the quantum cellular automaton A = {E,q, N,6), 
we denote by Ck...i the word Ck ■ ■ ■ ■ ■ ci if k < l, and the empty word e oth- 
erwise. Thus in either case Ck...i G (qU)* . Moreover we denote by Cí-^n the 
word Ci+^iniN)...i+ma^{N), and c^^ff the word Cí+inin(Ar)...í+max(íV)-i· Thcrcforc 
we have Cí+n G (<?-£')" a^f^ ^í+n ^ j respectively. 

Whilst configurations hold the bàsic states of an entire row of celis, and hence 
denote the possible bàsic states of the entire LQCA, the global state of a LQCA 
may well turn out to be a superposition of these. The foUowing definition works 
because C/ is a countably infinite set. 

Definition 3 (superpositions of finite configurations) 

A superposition of configurations of the quantum cellular automaton A = {S , q, N, d) 
is a normalized element ofTíCj , the Hilbert space of configurations. 

Definition 4 (global evolution) 

The global evolution of the quantum cellular automaton A = {E,q, N,S) is 
the linear operation defined by linear extension of its action upon the canonical 
orthonormal basis, as follows: 

A : Hcf — > l·lcf 
\c) ^ A\c) 

^|c) =(g)í|c,+iv) 

The postulates of quantum theory impose that the global evolution should be 
unitary. 

Definition 5 (Unitarity) 

The global evolution A of the quantum cellular automaton A — {S , q, N, 5) is 
said to be (finite) unitary if and only if {A\c) \ c G C/} is an orthonormal basis 
ofUcr 



The next three lemmas are known facts in linear àlgebra which will prové useful 
later on. 



Lemma 1 (Norm-preservedness, norm of rows) 

Let A : Tís — > Ti-s be a linear operator over a Hilbert space Tís having canonical 

orthonormal basis {\c)}ceS- Suppose the following two conditions are fulfilled si- 

multaneously: 

-(t)yc€S, [||Z\|c)|| = 1]; 

-(li) Vc, c' G S, [{c'\A''A\c) ^ ^ c = c']; 

ThenWeS, [0 < \\A^r)\\ < 1]. 



Proof: Conditions (i) and (ii) cxprcss the fact that A is norni-prcscrving. As a 
consequence for all r G 5 wc havc ||Z\'l'|r)|| = ||Z\Z\^|r)||. Moreover we have by 
definition ||Z\'l'|r)|p = | (r|/lZ\t |r) |. But the latter {{rlAA^r)) is a projection of 
the former (Z\Z\'l'|r)) over a unit vector (|r)), and hence |(r|Z\Z\'''|r)| < ||Z\Zi'''|r)||. 
Therefore ||^1'|r)||2 < \\A''\r)\\ and so \\A'i\r)\\ < 1. □ 
As a corollary we have: 

Lemma 2 (Norm-preservedness, finite spaces) 

Let A : Tís — > Ti-s be a linear operator over a Hilbert space Tís having canonical 

orthonormal basis {\c)}c£s- Let T be a finite subset of S. Suppose the following 

three conditions are fulfilled simultaneously: 

-(i)yceS, [||Zi|c)|r=l]; 

-(ii) Vc, c' G S, [{c'\A''A\c) ^ <^ c = c']; 

-(m)[EreT\\^Mr)\\^ = \T\]- 

Then we have Vr G T, [||Z\+|r)|| = 1]. 

Lemma 3 (Unitarity from unit rows) 

Let A : Tís Ti-s be a linear operator over a Hilbert space Tís having canonical 

orthonormal basis {\c)}ces- Suppose the following three conditions are fulfilled 

simultaneously: 

-r*;VcG5, [\\A\c)\\^í]; 

-(ii) yc,c' G S, [{c'\A''A\c) ^ <^ c = c']; 

-(ni) W G S, [\\A'<\r}\ \ = 1]. 

Then A is unitary. 

Proof: Conditions (i) and (ii) express the fact that {A\c) | c G S} is an orthonor- 
mal set, i.e. that A is norm-preserving. Condition (iii) expresses the fact that 
for all r G S, A^r) has unit norm. As a consequence Z\Z\^|r) has unit norm 
on the one hand, and (r|ziZi^|r) = ||zi^|r)|| = 1 on the other hand. Therefore 
Q;Z\Z\^|r) = \r), with a a root of unity. Since AA^ is positive, this a is just 1. 
Let A^\r) = X)/3c|c)- Then |r) = J2PcA\c), in other words each of the canon- 
ical orthonormal basis vectors |r) may be expressed as a linear combination of 
columns {A\c)\cG S}. Therefore the columns form themselves an orthonormal 
basis. □ 
Next we will examine each of the following conditions in turn: 



-(i) thc columns of A havc unit norm (section|3I); 
-(ii) the columns of A arc orthogonal (section^J; 
-(iii) thc rows of A have unit norm (scctionlS)). 



3 Unit columns 

The next two lemmas are simple facts from W. 

Lemma 4 (Norm of a column) Let A denote the global cvolution of the quan- 
tum cellular automaton A = [S, q, N, S). We have that 

WceCf, [||Zi|c)||=niHc,:+^)||] 

Proof. The norm of a tensor product of vectors is the product of the norms of 
the vectors. □ 

Lemma 5 ( Expressiveness of normalized S) Let A denote the global evo- 
lution of the quantum cellular automaton A = {S,q,N,5). Let A' denote the 
global evolution of the quantum cellular automaton A! = {S,q,N,S'), with S' 
such that Vw G (qS)", [S'\w) = (5|w)/| |(5|ïü) 1 1] . Suppose that the columns of A 
have unit norm. Then we have A ^ A' . 

Proof. For all c in C/ we have 

^|c> ^<^S\c,+n) ^ <^\\Sh+N)\\.S'\c,+N) 

= 11\\S\c,+n)\\.(^S'\c.+n) 

iei íGZ 

= \\A\c)\\.^5'\c,+r,)=A'\c) □ 

Our approach is to change the actual definition of LQCA as a consequence. 
Working definition 2 (LQCA) 

A linear quantum cellular automaton (LQCA) is a ^-tuple A = {S,q, N^S), 
where: 

- S is a finite set of symbols ("the alphabet"); 

- q is a symbol such that q ^ S ("the quiescent symbol"); 

- N is a set of n successive signed integers ("the neighboorhood" ) ; 

- S : Ti.(^qS}" ^ 'HqE ïs a function from superpositions of n symbols words to 
superpositions of one symbol words ("the local transition function"). 

Moreover S must verify the following two properties: 

- the quiescent stability condition: [S\q")) = jg)]. 

- the normalization condition: 
yweiqsr, [\\S\w)\\^l]. 



This modified definition of LQCA choses to imposo that Vw E (gZ")", [||<5|w)|| = 
l] , i.e. that the local transition function S is normalized. Then the fact that the 
columns of A have unit norm foUows straight from lemma^ We are strongiy 
justified to place this algebraic, almost syntactic restriction straight into the def- 
inition of LQCA because: 

- The alternative is to have various non-normalized states compensating each 
other non-locally, which from a physical point of view is somewhat disturbing; 

- It saves us from having to employ more elaborate techniques to check that 
columns have unit norms0|Zj, such as applying least path algorithm to the as- 
sociated de Bruijn graphs of the quantum cellular automata etc. Although very 
elegant these tend to render quantum cellular automata much more oblivious as 
a model of computation; 

- The modification made has absolutely no cost in terms expressiveness, as 
demonstrated in Lcmma|Sl 

4 Orthogonality of columns 

Having checked that the columns of the global evolution matrix A have unit 
norm, we now turn to the problem of deciding whether these columns are mu- 
tually orhogonal. First we need a definition. 

Definition 6 (A-matrix) Consider a linear quantum cellular automaton A = 
{E,q,N,S). We call A = [A'^y] with x,y e (gZ")""^ and a G (qS), the matrix 
(tensor) such that A^y = J2a ^xyW) equals S\w) if we have both x — íüi...,i_i 
and y = W2...n for some w e (qS)", otherwise it is the null vector. 

For instance consider the sample rule: (í|000) = |0), (5|001> = |1>, (5|010) = 
5\Q11) = |0), í|100) = |0), 5\1Q1) = í|110) = <5|111> = |0). 
Then the A-matrix of the sample rule is: 

/|0) |0) 0\ 

|1) |1) 

|1) |i) • 
V |0) 10)/ 

E.g. since 01 and 00 do not "foUow each other" the entry (00|A|01) holds the 
null vector. On the other hand 01 and 10 do overlap correctly to form the neigh- 
bourhood 010, for which í|010) = |1), hence (10|A|01) = |1). 
How can we check that all columns of some global evolution are mutually or- 
thogonal, when there is infinitely many of thcm? Our ncxt proposition is crucial 
in that respect, as it shows why the problem which might seem to be of an 
infinite (undecidable) nature is indeed of a finite (decidable) nature. 

Proposition 1 (Finite columns checks) 

This result refers to working definition [3 

Consider the global evolution A of a quantum cellular automaton A — {S , q, N, S). 
Let s be equal to jçi^p""^ — 1 and I be the interval [0, s]. The columns {A\c) \ c G 
Cf} are orthogonal if and only if the columns {A\c) \ c G Cj} are orthogonal. 



(See appendix lA. l·l for a detailed proof.) 

We now give our algebraic condition upon S ensuring that thc columns of A are 
mutually orthogonal. 

Proposition 2 (Column test) 

This result refers to working definition [21 

Consider the global evolution A of a quantum cellular automaton A = {S, q, N, S). 
Let s — Iq'Z'p"^^ — 1. The columns {Z\|c)|cgC/} are orthogonal if and only if 
yx,x' e {qS)"-\- 

[(a;a;'|M"|g"-ig""i)(g"-ig"-i|Ar|a;a;') ^ ^ {x ^ x')] (1) 

wüh M = [M^x'.yy'],M^^,^yy, = I T,a^x*y'A%y\^^ ^ the A-tensor of the LQCA. 
(See appendix lA. 21 for a detailed proof.) 

We believe that the obtention of algebraic conditions constituted a necessary 
step in ordcr to be able to take further the analysis of this model. The alge- 
braic proofs of these conditions give them a physical meaning which shortcuts 
the graph-theoretical detour, which is particular ly useful since quantum theory 
people tend to reason in terms of linear àlgebra rather than graph theory. They 
also master the corresponding numerical tools better, Proposition |21 makes it 
easy to check for column orthonormality through any software tool which does 
matrix multiplication. 

Still there remains some space for improvement, for instance because Proposi- 
tion 121 is phrased in terms of thc somewhat bizarre A-tensor of A, rather than 
just S. Fortunately this first point can be fixed using the quantum equivalent of 
a simplifying classical result [HjPT]jS]: 

Lemma 6 (Expressivity of size two 6) 

This result refers to working definition [21 

Consider a linear quantum cellular automaton A' = {S, q, N, S') and its global 
evolution A' . One can always construct a linear quantum cellular automaton 
A = g"^^, {0,1}, (5) such that its global evolution A equals A' . Then the 

A-tensor of A is just S — [S^y] , with x,y,a £ E"·~^ . 

Proof. We let, for aU x,y e Z""^^ 

6\xi . . . Xn-iyi ■ ■ ■ yn-i) = S'\xi...Xn-iyi---yi). 

The rest foUows from the definitions. □ 
Note that the groups of n — 1 celis constructed in this lemma bear some resem- 
blance with the reduced neighbourhoods N constructed by Dürr et al. - except 
they do not overlap, which consequently saves us from using De Buijn graphs 
and the like. This suggests that in spite of its apparent simplicity this trick is 
probably just the right way to enumerate/construct LQCA. 
We change the definition of LQCA as a consequence. 



Working definition 3 (LQCA) 

A linear quantum cel·lular automaton (LQCA) is a 3-tuple A — {S,q,S), where: 

- S is a finite set of symbols ("the alphabet"); 

- q is a symbol such that q ^ S ("the quiescent symbol"); 

- S : 'H(qS)'^ — > Ti-qE is a function from superpositions of 2 symbols words to 
superpositions of one symbol words ("the local transition function" ) . 

Moreover 6 must verify the following two properties: 

- the quiescent stability condition: \5\qq)) — \q)\. 

- the normalization condition: 
WweiqSr, [\\S\w)\\^l]. 

This modified definition of LQCA clioses to impose neighbourhoods of size two 
on top of tlie normalization condition. For these linear quantum cellular au- 
tomata we define A as usual with N = {0, 1}. Again we are strongly justified to 
place this easy restriction straight into the definition of LQCA for it simplifies 
and makes more intuitive the decision procedure induced by Proposition[21 (the 
corresponding simplifications are shown in CoroUary^. The modification made 
comes at absolutely no cost in terms expressi veness, as demonstrated in Lemma 

El 

These successive two modifications we have made to our model will be even more 
asserted by the important simplification they bring to the problcm of determin- 
ing whether a LQCA has unit rows. 

5 Unit rows 

Having chccked that the columns of the global evolution matrix A are orthonor- 
mal, we now turn to the problem of deciding whether its rows have unit norm. 

Proposition 3 (row norm as matrix product) 

This result refers to working definition |3| 

Consider a quantum cellular automaton A ~ {U, q, S) whose global evolution A 
has orthonormal columns. The squared norm of any row r is given by 

||ZiV)||2^ lim (g|iV('?)''ï TT Ar('-'))iV(9'''|g) 

iek...i 

where {N^''')}^çqE is the set of matrices such that N'^''^ = [ivifj]], Tvifj) = \{a\S\xy)\'^ . 
(See appendix lA.3l for a detailed proof.) 

The following proposition takes advantage of the successive restrictions made in 
definitionsElandinito bring about a crucial simplification - which makes obsolete 
a good half of the procedure described in 



Proposition 4 (Middle segment) 
This result refers to working definition|3| 

Consider a quantum cel·lular automaton A = {U, q, S) whose global evolution A 



with O = [Ixy] the matrix with only ones, and N = [Nx,y], N^ y — \{q\5\xy)\'^ . 
(See aDT)endix lA.4l for a detailed proof.) 

Note that \iinh^^{q\N''ON''\q) = V ^^[o,oi ||Z\'''|r)|P. Hence we have the fol- 

lowing insightful coroUary, which comes as the direct analogue of our Proposition 
n Curiously however it will not contribute to our final result. 

Proposition 5 (Finite rows check) 

This result refers to working definition O 

Consider a quantum cellular automaton A — {U, q, S) whose global evolution A 
has orthonormal columns. The rows {^^|?') | r S C/} have unit norm if and only 
if the rows {A^r) \ r £ Cy*'"'} have unit norm. 

Both results deepen our understanding of the algebraic structure of unitary lin- 
ear quantum cellular automata. Proposition ^offers an important simplification 
along the way to determining whether the global evolution A has unit rows, 
reducing this to the evaluation of l .0 r . 

The difficult problem we are left with is that of evaluating the so-called 'border 
vectors' [5] "r* = lim?i^oo N'^\q) and l = hm/j^oo TV''' \q). The foUowing propo- 
sition will characterize them uniquely and algebraically. We let < denotes the 
foUowing partial order upon m x n matrices: 

M < N yi,j < N,j. 

Column vectors are seen as m x 1 matrices for that matter. 

Proposition 6 (Border vectors) 

This result refers to working definition|31 

Consider a quantum cellular automaton A — {S, q, 5) whose global evolution A 
has orthonormal columns. The vectors 



with N = [Nx.y], Nx.y = \ {q\5\xy)\'^ , have only finite entries. They verify 




r = 




lim N''\q) and T = lim N^''\q) 



= minív I < u A Nv 

< 

and l = minju \ < v A vN 



V /\ Vq = 1} 



V /\ Vq ^ 1}. 



Moreover the foUowing extra conditions hold: 
(i) t.V = 1; 



(Ui) Va; e [T^ = V T*^ = 0]; 

(w) Vx, yeS, [N,y ^ T, = V Ty = 0]; 

Inequality (ii) is saturated if and only if A has unit rows. 

(See appendix I A . 61 for a detailed proof.) This last proposition is highly informa- 
tive, and in most cases will provide us with an efFective way to compute these 
border vectors, through a spectral decomposition of N. When the eigenvalue 1 
is degenerate, however, it is unclear to the author whether there exists a definite 
procedure to performing the minimization. 

Taking a step back from the mathematics, we may wonder where these hmits 
N'^, /i — > oo come from 'physically'. Say we wish to calculate the norm of a 
row Z\^|r), corresponding to a configuration r with interval domain /. Then at 
some point we need to sum over all the antecedent configurations of r (cf. the 
beginning of the proof of Proposition|31). Our issue now arises from the fact that 
although the antecedents of r are in C/, there may be an infinite number of them, 
and hence they may not be able to restrict to some interval J. A good example 
of this is provided in which we now reproduce for convenience: 

Remark 1 Consider Qflip {p,q,S) with S\qq) = \q), S\qp) = \q) + |p)/·\/2, 
S\pq) = \p), 6\pp) — \q) ~ \p)/V2. Its corresponding global evolution A is unitary. 
Let c" — . . . qqp"qq . . ., with the rightmost p in position 0. Then (c^|Z\|c") = 
(1/a/2)". 

Faced with such situations our temptation is again to further restrict the defi- 
nition of LQCA, with a view to eliminate such scenarios. For instance we may 
want to add to the Working definition |31 a full stahility condition asking that: 
VwG í:^^ [{q\S\w) =0]. 

This extra 'full stability condition' would ensure that the antecedents c of a 
configuration r have an interval domain no greater that the extended interval 
domain of r. 

Lemma 7 (Full stability unit rows) 

This result refers to working definition |2| 

Consider a quantum cellular automaton A — {S, q, 5) whose global evolution A 
has orthonormal columns. Suppose the full stability condition\/w G S'^ ^ \{q\5\w) = 
O] . is verified. Then the rows {A^r) \ r G Cb} have unit norm if and only if 

{q\NON\q) = \qE\ (2) 

with O — [Ixy] the matrix with only ones, and N — [Nx^y], Nx^y = \ {q\S\xy)\'^ . 



Proof. First note that undcr this condition, one has "ix £ S 



{x\N^\q)= N.,.N^,^ Yl \mx<y)\'\{q\S\aq)\' 

= \{q\S\xq)f\(q\S\qq)f by the full stability condition 
= |((7|(5|a;g)p by the quiescent stability condition 
= {x\N\q). 

Moreover we have {q\N^\q) — J2aeqS -^qi^^crQ^ which is equal to NqqNqq = 1 = 
{q\N\q) using the fact that Vx G i7 [NqxNxq = 0]. The contrary would imply 
{. . . qq . . . \ A\ . . . qqxqq . . .) 7^ 0, which is impossible since {q\5\qq) = 1 and A is 
norm-preserving. 

As a consequence we have N'^\q) = N\q). Symmetrically so for (çliV^ = {q\^- 
The rest foUows from Proposition^] □ 
Unfortunately it is not clear whethcr this extra restriction of LQCA implies a 
loss of expressiveness. In fact may even disqualifies some vàlid reversible cellular 
automata from being quantum cellular automata. We leave the question of a 
simulation as an open problem. 

Open problem 1 Consider a linear quantum cellular automaton A — {S, q, S) 
following to definition\^ and its global evolution A. Is it always possible to con- 
struct a linear quantum cellular automaton A! = (Z", g', 5') such that Vw € 
[((7|(5|w) = 0] and its global evolution A' equals A? (We may want to admit 
loser notions of ejjicient simulation.) 

In fact one may wonder whether there may not be a way to make any norm- 
preserving LQCA (unit, orthogonal colunins) into a unitary LQCA exhibiting 
the full stability condition. For instance: 

Remark 2 Consider XOR = ({0,1}, 0,(5) with 5\00) = |0), S\01) = S\10) = 
\1) , 5\íí) — |0). Its corresponding global evolution A is injective in the space 
of finite configurations (i.e. it is norm-preserving, it has orthonormal columns). 
But the rule S\ll) = |0) breaks the full stability condition. Moreover it is not 
surjective (i.e. its rows are not unit norm) since . . .00100. . . does not have any 
antecedent. 

is suitably 'fixed' as in Remark 01 
6 Conclusion 

The following theorem is a synthesis our results. 

Theorem 1 (Summary) 

This result refers to working definition |2| 

Consider the global evolution A of a quantum cellular automaton A = {S, q, S). 



Let s — \qS\'^ — 1. A is unitary if and only if\fx,x' £ (q^)-' 



[{xx'\M''\qq){qq\M'\xx') ^ ^ {x ^ x')] 
and l i)ÇYl^i) = \qS\ 

i i 

with ~r* = min{í; | < v /\ Nv = v A Vq = 1}; 
l — min{í; | < w A vN — v A Vq — 1}; 

N^[N^,y], N^,y^\{q\5\xy)\^. 

We have definitely gone a long way towards the simplification and the alge- 
braization of unitarity criteria for one dimensional quantum cellular automata 
as defined in [S]. Note that these last two papers do not contain any 

such synthetic algebraic criteria. Instead they provide several pagès long decision 
procedures, which have many twists and bends. 

Certainly we have not gone as far as to reduce LQCA to partitioned quan- 
tum cellular automata (PQCA) or quantum cellular automata with Margolus 
neighbourhood (MQCA) etc, i.e. our global evolution cannot, in general, be 
decomposed into the application of one small unitary operator homogeneously 
across space. 

Other models of quantum cellular automata do admit such reductions; the prob- 
lem of deciding unitarity becomes incomparably easier, or even trivial by con- 
struction ^12, . Such crucial differences seem to arise when one considers different 
spaces of configurations, e.g. finite periòdic because these constrain re- 

versibility to be structural, i.e. an essentially local matter. With the space of 
finite configurations Cf reversibility becomes a global matter, even though the 
global evolution is defined locally. Anàlogues of this are well-known in the clas- 
sical realm already: 

Remark 3 Consider XOR' = ({g, 0. 1}, ç, (5) with 5\00) = |0), S\OÍ) = \í), 
S\ÍO) = S\n) = |0),, 6\qO) = \q), S\ql) = \q) , S\Oq) = |0), S\lq) = 

Its corresponding global evolution A is bijective in the space of finite configura- 
tions (i.e. it is unitary), but the global evolution cannot be reversed by a cellular 
automata. Moreover A is not reversible in the space of infinite configurations 
(. . .00 . . . may have antecedents either . . . 00 . . . or ... 11 .. .j. 

One may argue that PQCA models are enough, since they simulate the quan- 
tum Turing machine. If we are interested in 'intrinsic universality' however, and 
want to consider the above example as 'physical' then these are not be enough. 
A local transition S which is unitary in the space of infinite configurations is also 
unitary in the space if finite configurations, but the reciprocal statement is not 
true. In this sense there seems to be a loss of expressiveness when restricting to 
PQCA models. 



Note that this research constitutes a first step in the quest for the identifica- 
tion of the onc particular unitary LQCA known to simulatc all othcr LQCA 
efficicntly. Obviously as we restrict the number of one dimensional quantum cel- 
lular automata to be considered to a well-behaved subclass, the quest for such 
an 'intrinsically universal' onc dimensional quantum ccUular automata must be- 
come easier. Moreover our approach may eventually opcn this (well-advcrtised 
in the literature) opcn problem to an algcbraic analysis. 

Note also that it is undecidable whethcr local transition function S induccs a 
unitary global evolution A as soon as the quantum cellular automata is two 
dimensional. But again this does not havc to bc the end of it: nothing prevcnts 
that thcrc should bc a canonical definition of two dimensional quantum cellular 
automata, easily chccked for unitarity, and yct capablc of cxprcssing the cxact 
same set of global evolutions. For this purpose algcbraic criteria ought to be 
easier to generalize to higher dimensions. 
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A Proofs of Propositíons HHül 
A.l Proof of Proposition ^ 

Consider c 7^ c' such that A\c) is not orthogonal to Zi|c'). Since there are Igi^p""^ 
possible pairs of (çZ")""^ subwords, there must exist k,l ^ Z with < l — k < s 
such that Ck...k+n-2 = Q...i+n-2 and 4 j,^„_2 = 4...i+n-2- We now construct 
d = Coo...k+n-2 ■ Ci+„-i...oo and d' = c'^,„k+n-2 ■ cí+„_i...oo' and notice that 

{A\d),A\d'))= n {5\c!,+^).5\c,+m))^0 

í^]fe+?i-2,í+ri-l[ 

Since {A\c),A\c')) = [] {5\4+j,),5\c,+n)) ^ 0. 

i 

By iterating this procedure the right number of times we obtain, from any c 7^ c' 
such that Z\(c) is not orthogonal to ^(c'), a pair d ^ d' such that A{d) is not 
orthogonal to A{d') and with both idom((i) and \àom{d') included in an interval 
of size less or equal to s. By shift-invariance it is therefore sufíicient to check 
that the columns {A\c) | c € A idom(c) < s} are pairwise orthogonal. □ 



A.2 Proof of Proposotion El 



Let / stand for the interval [0, s] and J stand for the interval [— s, s]. 
[=>] CoroUarynimplies that the columns {Z\|c) |c e C/} are orthogonal if and 
only if we have Vc, c' G | (Z\|c), zí\|c')) p strictly positive when c equals c' 
and zero otherwise. : 



c—dl c^d 

iGcxtidom( J) 
'íGcxtidom( J) 
zÇextidom( J) 

^ TT - c' > / = 

2í^cxtidom( J) ^ 



Summing this last equation over those configurations c, c' G which verify 



^o+JV = ^' "^n^w ~ ^' yields, for any fixed x, x' € 17" ^: 



{q^'-^q^-'^\M'\xx'){xx'\M'\q''-^q''-^) > /=0 

x=x' I x^x' 

Consider two configurations d, d' G C^. By merely shifting to the center any 
difference between d and d! one can always construct c, c' € Cj such that Cq_|_^ = 
X, c'qj^j^ — x' , and [x = x' if and only if d = d']. Now consider condition 



■i 

■Mxx',...\\M, . -,c.^,^^c' _M...,ç„-i,„-i > /=0 

11 ^ + N i + jv' ! + l + N i+i + JV ''^ X=X'IX^X' 

i 

^ n m+N))\5\c[^^)) >J=Q 

c—dl C^C 

iGextidomí J) 

^ n ('5M.+iv))n^M::+^» /= 

-^-^ d=d'/ d^d' 

zGcxtidoin(/) 



In other words the columns {^|<i) | d € Cj} are orthogonal, and so CoroUary^ 
imphes that the columns {A\c) \ c S C/} are orthogonal. □ 



A.3 Proof of Proposition |31 



r,ceCf 

= E \r){r\^S\c,c,+i){c\ 
r.ceCf ieï 

ceCf iez 

{r\AA'^\r)=Y^\{\{n\5\c^c,+i)\'' 

ceCf iez 

ceCf iez 

Say r has an interval domain included in / = [k,l]. Thcn ||Z\^|r)|p equals 

E n m,.. n n 

cÇC/ zíE— 00.../C — 1 i{^k...l 

If we now restrict our sum to antecedents whose interval domain is included in 
J — [k — h,l + h + l], cach tcrm in the sum takes the form 



gq qCk-h Cfc_iCfc V J_J_ CiCi^il Ci^iCi^^ Ci^^^iq qq 

iek...i 

= N^i) . . . mM ( Y\ Min) ^/Y(«) ...Mi'l) 

i£k...l 

using Nj^^ — 1. Performing the sum yields 

fc-l-/i...fc-l iek...l l+l...l+h+l 



■iGfc...i 

The limit of this expression as h tends to infinity remains finite. Indeed, let us 
suppose that the contrary is true. This implies that the row A^r) has infinite 
norm. But by lemmaHthis contradicts the fact that A is assumed to be norm- 
preserving. □ 

A.4 Proof of Proposition |4] 

By lemma|21 requiring that the rows be of unit norm of rows is equivalent to the 
condition that for all t and J = [— í,í], 



E \\A'^\r)\? = \qS\ 



2Í+1 



rec/ 



Using proposition|31thc cquation becomes 



lim {q\N''{ ^ '[[N^''^^)N''\q) = 

recj ieJ 

^ lim {q\N'^{ Y[ ^ 7V(-))iV"|g) = \qS\''+' 

-Í...Í o-eqU 

^ lim {q\N''0^'+^N''\q) = \qS\'^'+^ 

h — >oo 

since ^ Nil^= ^ = 1 

by the normalization condition. 
Since O^t+i _ have our condition. □ 

A.5 Proof of Proposition |5] 

Since \imii^oo{q\N'^ON'^\q) = X^rec'""' H'^^l'')!!^ ^^^^ comes as a coroUary of 
Proposition 21 D 



A.6 Proof of Proposition |ül 

— )■ ^ — > — > 

We have that l . r — l q — l q — 1 since all three are by definition equal to 
lim/i^oo(<z|^''|'í) = ll^^l ■ • • 99<Z ■ • ·)IP = l-ÍThat tlic norm of thc all-quiescent 
row is one stems from the fact that {q\S\qq) ~ 1 and A is norm-preserving.) 
All entries in l and are nonnegative because all entries in N and \q) are 
nonnegative. 
Now note that 

(E'^^)(E^'0= lim {q\N'^ON'^\q) 
i i 

= J2 \\^Mr)\\'<\qS\ 

where the last line was obtained using J^aeqS ^^'v ~ ^ ^^'^ Lemma (It 
is clear that this inequality is saturated if and only if A has unit rows, since 
this was the content of proposition 01) Using this and nonnegativity we have 

(J2i^i) = ^ qiY^i^í) < (I2i l i)iY.i^i) < k^l' lience the entries in are 
finite, and symmetrically so for l . 

Saj 3x G S, [tx = /\ = 0]. Then {. . . qq . . .\A\ . . . qqxqq . . .) ^ 0, 

which contradicts the fact that {q\S\qq) = 1 and A is norm-preserving, and 

yields condition (iii). Say 3x,y e S, [N^y ^ A l x — ^ ^^==0]- Then 
{. . . qq . . . \ A\ . . . qqxyqq . . .) 7^ 0, which again contradicts the fact that {q\5\qq) = 
1 and A is norm-preserving, and yields condition (iv). 



The fact that is an eigenvector of N is trivial, since = lim;i^oo -^''k) = 
lim/i^oo N'''^^\q) — NT^. We now show it is minimal in 5 = | < i; A Nv = 
V A Vq = 1}. Say w Cz S. We have w = \q) +u, with u a vector liaving nonnegative 
entries. Hence 

w = lim N'^w 

h — ^oo 

= lim N''\q) + lim N'^u 

h^oo h — >-oo 

= y + lim N'^u. 

h — >oo 

But sincc N and u have nonnegative entries, lim/i^oo -/V^m has nonnegative en- 
tries and ~r' < w. Symmetrically so for l . □ 
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